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VOLUME GROWTH AND ESCAPE RATE OF BROWNIAN 
MOTION ON A COMPLETE RIEMANNIAN MANIFOLD 

By Elton P. Hsu and Guangnan Qin 

Northwestern University and Chinese Academy of Sciences 

We give an effective upper escape rate function for Brownian 
motion on a complete Riemannian manifold in terms of the volume 
growth of the manifold. An important step in the work is estimating 
the small tail probability of the crossing time between two concen- 
tric geodesic spheres by reflecting Brownian motions on the larger 
geodesic ball. 

1. Introduction. Let M be a Riemannian manifold and pM(t,x,y) the 
(minimal) heat kernel on M. By definition, the latter is the fundamental 
solution of the heat operator 

^ M = Wt~\^ Ah 

where Am is the Laplace-Beltrami operator on M. A Riemannian manifold 
is stochastically complete if 

p M (t,x,y)dy = 1 

M 

for some, hence all, (x,t) S M x (0,oo). In other words, M is stochastically 
complete if the heat kernel is conserved. Let F x be the law of Brownian 
motion on M starting from x and let e be the lifetime (or explosion time) 
of the Brownian motion. We then have 

F x {e>t}= / p(t,x,y)dy. 
JM 

Therefore, M is stochastically complete if and only if 
(1.1) P a .{e = oo} = l, 
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that is, Brownian motion on M does not explode. Finding geometric con- 
ditions for stochastic completeness is an old geometric problem. The prob- 
lem has been attacked using both analytic and probabilistic methods. Early 
works on this problem (Karp and Li [9], Yau [16], Hsu [6] and Varopoulos 
[15]) impose lower bounds on the Ricci curvature. In particular, in the last 
two works, it was shown that if there is a strictly positive function n(r) such 
that the Ricci curvature of M on the geodesic ball B(r) is bounded from 
below by — «(r) and 



then M is stochastically complete. In 1986, Grigor'yan [1] found the following 
sufficient condition for stochastic completeness solely in terms of the volume 
growth function of the manifold: 



According to (1.1), a Riemannian manifold M is stochastically complete 
if Brownian motion does not escape to infinity (in the one-point compactifi- 
cation) in finite time. Traditionally in probability theory (see, e.g., ltd and 
McKean [8] and Shiga and Watanabe [12]), one often looks for upper func- 
tions for the escape rate of a diffusion process. The classical Khinchin law 
of iterated logarithm for one-dimensional standard Brownian motion is the 
most celebrated case. More generally, let r t = d(X t ,x) be the radial process 
of Brownian motion X on M. An increasing function ip(t) is called an upper 
rate function if 



Finding an upper rate function is a more refined problem than proving 
stochastic completeness since the existence of an upper rate function implies 
stochastic completeness. Various explicit upper rate function for Brownian 
motion on a complete Riemannian manifold have been obtained under con- 
crete volume growth assumptions (see Grigor'yan [2], Grigor'yan and Kel- 
bert [4] and Takeda [13, 14]). More recently, Grigor'yan and Hsu [3] showed 
that the inverse function of the increasing function 



related to the integral test (1.2) for stochastic completeness is essentially 
an upper rate function for Brownian motion on M. While this result gives 
an upper rate function of a very general form, it was proven under the 
additional geometric assumption that M is a Cartan-Hadamard manifold, 




(1.2) 




^x{ r t < i>(t) for all sufficiently large t} = 1. 
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that is, a simply connected, geodesically complete Riemannian manifold with 
nonpositive sectional curvature. 

The main purpose of the present work is to obtain an escape rate function 
based solely on the volume growth of the underlying manifold. We introduce 
the following increasing function: 



We will show in Theorem 4.1 that under the sole assumption that M is a 
complete Riemannian manifold, the inverse function of (j) is essentially an 
upper rate function for Brownian motion on M. 

Remark 1.1. The difference between the functions cpi(r) and (f>(r) is 
that we have introduced an extra term lnlnr* in the latter function. This 
addition, resulting from an attempt to remove the extraneous geometric 
condition, is fully justified on several grounds. First, the integral test (1.2) for 
stochastic completeness, which can be written as ^i(oo) = oo, is equivalent 
to the condition 4>(oo) = oo. Second, our new upper rate function implies all 
explicit upper rate functions existent in the literature to date (see Corollaries 
4.2 and 4.3). Third, in general, the radial process of a Brownian motion on 
M has the form 



where (3 is a standard one-dimensional Brownian motion and L is a local 
time on the cut locus C(x) of the point x (see Hsu [7] and Kendall [10]). In 
the absence of any further geometric assumptions, we do not expect to obtain 
an upper rate function (up to a multiplicative constant) for the process r t 
better than the upper rate function of a standard Brownian motion ijj(t) = 
Cy/i In Int. This rate function cannot be achieved without the presence of 
the additional term lnlnr* in the function <j>. 

Our method has two key steps. In the first, we follow Hsu [6] and Grigor'yan 
and Hsu [3] and, by using the Borel-Cantelli lemma, reduce the problem of 
seeking an upper rate function to the problem of estimating the small tail 
probability of the crossing time between two concentric geodesic spheres 
(Lemma 2.1). In the second key step, instead of estimating the small tail 
probability by using an analytic approach, as in Grigor'yan and Hsu [3], 
under the assumption that the manifold in question is Cartan-Hadamard, 
we modify the way the Lyons-Zheng decomposition for reflecting Brownian 
motion is used in Takeda [13, 14]. The volume |-B(r)| of the geodesic ball 
B(r) appears naturally in this step because the uniform distribution on the 
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ball with respect to the Riemannian volume measure is the invariant mea- 
sure of reflecting Brownian motion on B(r). The additional term mentioned 
above, lnlnr, is a consequence of dealing with the Brownian motion adapted 
to the time-reversed filtration in the Lyons-Zheng decomposition. 

2. Basic estimates on crossing times. Let M be a geodesically complete 
Riemannian manifold and &(M) the path space over M. Let X be the 
canonical coordinate process on the path space ^(M) over M, that is, 
Xt(u) = u>t for oj E ^(M). If x S M, then we use ¥ x to denote the law of 
Brownian motion on M starting from x. The radial process is rt = d(Xt,x), 
that is, the Riemannian distance from x to Xt, the position of Brownian 
motion at time t. A nonnegative increasing function R:M. + — > K + is called 
an upper rate function for Brownian motion on M if 

^x{ft < R{t) for all sufficiently large t} = 1. 

Let {R n } be a strictly increasing sequence of positive numbers, to be 
chosen later, and define a sequence of stopping times as follows: 

T n = ml{t :r t = R n }. 

(We adopt the convention that inf = oo.) Each of these is the first time 
the Brownian motion X reaches the corresponding geodesic sphere, 

S(R n ) = {x£M: d(x, o) = R n }. 

The difference r n — r n _i, if well defined, is the amount of time Brownian 
motion takes to cross from S(R n -i) to S(R n ). The basic idea of Grigor'yan 
and Hsu [3] (see also Hsu [6]) for controlling the rate of escape of Brownian 
motion is to give a good upper bound for the small tail probability T x {r n — 
T n -i < t n } for a suitably chosen sequence {t n } of time steps. If the sum of 
these probabilities converges, then the Borel-Cantelli lemma shows that for 
sufficiently large n, Brownian motion X has to wait roughly until at least 

n 

T n = tu 

k=l 

to reach the sphere S(R n ) or, equivalently, rt < R n for all t < T n . This, after 
some technical manipulations (see Section 4), will give an upper escape rate 
function. 

We now use the idea of Takeda [13, 14] to estimate the small tail prob- 
ability P x {r n — r n _i < t n } by using the Lyons-Zheng decomposition [11] 
of reflecting Brownian motion starting from the uniform distribution on a 
geodesic ball. For an open set B C M, we denote by Fb the law of Brownian 
motion starting from the uniform distribution on B, that is, 

r B = — I p x dx. 

\B\Jb 
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Likewise, we use to denote the law of reflecting Brownian motion on 
B starting from the same uniform distribution. Let B n = B(R n ) be the 
geodesic ball of radius R n centered at x. In order to take advantage of the 
volume growth condition, we consider the probability F Bl {r n — r n-i < tn} 
instead of ¥ x { T n ~ T n-i < t n }. Recall that r n is the first time the process X 
reaches the boundary S(R n ) of the geodesic ball B(R n ). Before reaching the 
boundary, Brownian motion and reflecting Brownian motion have the same 
law. Therefore, if C S SS Tn is an event which is measurable up to time r n , 
then ¥ Bn (C)=Q Bn (C). From 

f b 1 (C) = 7± 1 [ P x (C)dx<^- f F x (C)dx= l -^\¥ Bn (C), 

\ B l\ JBx \B\\ Jb„ \Bl\ 

we have 

(2.1) rBii c)<\^\Q Bn (c), CeJ T „. 

We apply this inequality to the event 

(2.2) C = {T n -T n _ 1 <t n }. 

Now, according to the Lyons-Zheng decomposition [11], on a fixed time 
horizon [0,T n ], the radial process can be decomposed as the difference 

Bf Bt — Bf 

(2-3) n -r = -±- T " 2 \ 

where B is a standard Brownian motion adapted to the natural filtration 
msi, = &(&(M))* of the path space &(M) and B is also a standard Brow- 
nian motion, but adapted to the reversed filtration ^* defined by 



a{X Tn _ s :0<s<t}, 0<t<T n . 



The advantage of such a decomposition is obvious, for we have eliminated 
from consideration the bounded variation component of the radial process, 
which can be rather complicated. The price is that we have to deal with a 
Brownian motion not adapted to the original filtration. Another complica- 
tion is that the decomposition cannot be applied directly to the event (2.2) 
because it may go beyond the fixed time horizon [0, T n ]. In order to remedy 
this situation, we will use a slightly modified event 

in the inequality (2.1). Fortunately, this additional restriction {r n < T n } will 
not be an obstacle for us, as shown by the following simple observation. 

Lemma 2.1. Let C n be defined as above. Suppose that X^nLi^(Cn) < °°- 
Then, with probability 1, there exists T_\ such that T n >T n — T—\ for all n. 
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Proof. By the Borel-Cantelli lemma, the probability that the events 
{C n } happen infinitely often is 0. Therefore, with probability 1, there exists 
no such that for all n > no, either r n — r n _i > t n or r n > T n . We show, by 
induction, that T n >T n — T no holds for all n. If 1 < n < no, then r n > > 
T n — T no . Suppose that r n > T n — Tq for an n > no- If r n+ i > T n+ i, then, 
trivially, r n+1 > T n+1 - T no . Otherwise, r n+ i - r n > t n+1 and 

Tn+l = T n+ l — T n + T n > t n+ \ +T n — T nQ = T n+ \ — T nQ . 

This completes the proof. □ 

We now prove the main estimate for the crossing time r n — r n _x- 

Proposition 2.2. Let r n be the first hitting time of the sphere S(R n ) 
and r n = R n — R n -\. There then exists a constant C such that 

F Bl {Tn - r„-i < t n ,r n < T n } < °^ e -^. 



\Bi\ y/irtn r n 

Proof. The event {r n — r n _i < t n } implies the event 
( 2 - 4 ) { sup (r Tn _ 1+s -r Tn _ 1 ) > r n \. 

Now, from the decomposition (2.3), we have 

2(r Tn _ 1+s - r Tn _ 1 ) = B Tn _ 1+s - B Tn _ x + B Tn _ 1+s - B Tn _ x . 

Since r n _i is a stopping time with respect to the natural filtration £8*, the 
first term on the right-hand side is a Brownian motion in time s starting 
from 0. This is not so for the second term because r n _i is not a stopping 
time with respect to the filtration of the reversed process. However, for 
any s < t n such that r n _i < T n , taking k such that (k — l)t n < r„_i < kt n , 
we see that both r n _i and r n _i + s lie in the interval \{k — l)t n , (k + l)t n ]. 
From 

r T„-i+s ~ r T„_i = r r„_i+s — r kt n + r kt n ~ r r„-i i 

the event (2.4) is contained in the union of the [T n /t n ] + 1 events 

X" 



sup \r ktn+s -r ktn \ > — }, l<k< 

\s\<t n 



+ 1. 



Using 
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we see that the event {r n — r n _i < t n , r n < T n } is also contained in the union 
of the following 2[T n /t n ] + 2 events: 



and 



sup \B ktn+s 

\s\<t„ 



sup \B ktn . 

\s\<t„ 



ft-i'n I — 2 



B kt n | > 



for 1 < k < [T n /t n ] + 1. Under the probability Qb„, these events have the 
same probability, 



mp\B tn+s -B tn \>^}<2l 



s\<t n 



sup \B S \ > — - } < 

0<s<t n ^ 



C\ft n 



The probability Q,B n {T n — r n __i < i n ,T n < T n } is bounded from above by 
2[T n /t n ] + 2 < 4T n /t n times the above probability. The desired inequality 
now follows immediately from this and the inequality [see (2.1)] 

\B n \ 



FBi{r n - r n _i < t n , r n < T n } < 



B, 



IB, 



,\j~n T"n—1 ^ t n ,T n < T ri \. 



□ 



3. Total crossing time. In the preceding section, we have found an upper 
bound for the probability Pb^t™ — r n _i < t n , r n < T n }. We are still free to 
choose the upper bounds t n of the crossing times r n — r n _i and the radii R n 
of the expanding geodesic balls B(R n ). We need to choose them so that the 
series 



oo 

£ 

n=l 



Fsi{r n - r n _i < t n , T n < T n } 



C 



E 

n=l 



^nXn 



exp 



hx\B(Rn)\ - 



r< 00 rp 

^LBd2^^ eXp 

111 n=l n 



ln\B(R n ) 



-2 "I 



16t„ 



converges and the Borel-Cantelli lemma can be applied. The obvious choice 
is for t n equal to a small multiple of r^/ln\B(R n )\ , as was adopted in 
Grigor'yan and Hsu [3]. However, this choice will not enable us to eliminate 
the extra factor T n jr\, whose presence can be traced back to the Brownian 
motion B adapted to the reverse filtration in the Lyons-Zheng decom- 
position (2.3). We diminish the obvious choice by letting 

1 r 2 

{ ' ] n Z2\n\B{R n )\+h{R n ) 
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with a strictly increasing function h to be determined. If we assume, without 
loss of generality, that B(R\) > 1 and h{R\) > 1, then t n < r^/32. If we 
further assume that the sequence {r n } is increasing, then there is an obvious 
bound 

n n 

32T n < r\ < r k r n = R n r n . 
k=i k=l 

It follows that 

P Bl {r n - r B _i < t n ,r n < T n } < ^^e~ 2h ^. 

\Bi\ 2 r n 

It remains to choose the radii R n and the function h such that 

OO OO 

(3.2) J>„ = oc and ^ f^e" 2 ^") < oo 

n=l n=l rn 

under the integral condition 

. . f°° rdr 

The divergence of the above integral is to be linked to the divergence of the 
total crossing time in (3.2). This leads to the natural requirement that 

r\ > CR n (Rn + i — R n ). 

This requirement can be fulfilled by setting R n = 2 n with C = 1/4. From 
(3.1), we have 

f o a) T - V — V R k( R k+i-Rk) 1 f Rn+1 rdr 

U ~ ~ 128 ^ ln|J3(i2fc)| + h(R k ) ~ 256 J Rl ln\ B (r)\ + h(r)> 

which seems to fall slightly short of the condition (3.3). The apparently 
disadvantageous situation can be salvaged by first looking at a typical can- 
didate for the function h. From our choice of R n — 2 n , we have R n j r n — 2 
and the convergence of the total probability in (3.2) becomes 

oo 

J2e~ 2h{2n) <oc. 

n=l 

This leads to the choice h(R) = lnlni?. We have the following simple obser- 
vation. 

Lemma 3.1. Let f be a positive, nondecreasing and continuous function 
on [0, +oo) such that 

(3 ' 5) J, 7M = °°' 



VOLUME GROWTH AND ESCAPE RATE OF BROWNIAN MOTION 



9 



Then 

rdr 



fir) + lnlnr 



■ oo. 



Proof. Divide the integral into the sum of the integrals over the inter- 
vals [re — l,re] for n > 4. Since / is increasing, we have 

rdr ^ n — 1 

f(r) + lnlnr ~~ fin) + In In re 

n=4 v ' 

1 x - re — 1 1 \ - n — 1 
- 2 ^ f(n) + 2 ^ In Inn" 

/(n)>lnlnn J v ' /(n)<lnlnn 

Since (re — 1)/ In In re > 1 for all sufficiently large n, if the second sum has 
infinitely many terms, then it is clearly diverges; otherwise, /(re) > In Inn 
for all sufficiently large n and we have, for some no, 



rdr 1 A re-1 



3 



It, 



f(r) + lnlnr 2^ /(n 

1 ^ n - 1 /- nH 
2" 4- reTll 



n=no 

rdr 



> 

2 

n=no 

1 n - 1 f 00 rdr 



> 



2 n + 1 J no f{r 



This completes the proof. □ 



We are now in a position to bound the range of Brownian motion on a 
finite time interval. 

Proposition 3.2. LetR n = 2 n and 

1 R k (R k+1 -R k ) 



1 98 ^ 



128f^ln\B(R k )\ + h(R k ) 

Then, with probability 1, there exists T—\ such that sup <£<r — t_i r * — 2 n 
/or a// n. 

Proof. By our choice of R n , 

oo 

/^B^n ~ Tn-1 < t n ,T n < T n } < OO. 
n=l 
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By the Borel-Cantelli lemma and Lemma 2.1, with probability 1, there exists 
T_i such that r n > T n — T-\. However, r n is the hitting time of the sphere 
S(2 n ), hence 

su Po<t<T n -T_i r t — 2 n for sufficiently large n. □ 

An easy consequence of the above result is a probabilistic proof of Grigor'yan's 
criterion for stochastic completeness. 

Corollary 3.3 (Grigor'yan [1]). Suppose that M is a complete Rie- 
mannian manifold and B(R) its geodesic ball of radius R centered at a fixed 
point. If 

rdr 

oo, 



ln|fl(r)| 

then M is stochastically complete. 

Proof. By Lemma 3.1, under the assumption, we have 

1 r Rn+1 rdr 
n ~ 256 J Rl ln\B(r)\+h(r) ~* °° 

as n — > oo. By the above proposition, sup t<T r t < oo for all finite T. Hence, 
Brownian motion does not explode and M is stochastic complete. □ 

4. Upper rate function. Proposition 3.2 allows us to obtain an upper 
rate function in terms of the volume growth function |U(r)|, as was similarly 
done in Grigor'yan and Hsu [3]. 

Let 

m fR rdr 



Yn\B{r)\ + lnlnr 

From (3.4), we have (l/256)(/>(2 n+1 ) <T n . Proposition 3.2 then gives 

sup r t < 2 n 

<<(l/256)</.(2"+ 1 )-T_ 1 

for all n > 1. This implies that 

sup r t < 2R 

i<(l/256)^(fl)-T_! 

for all R > 0. Denote by ip the unique inverse function of </>. Letting R = 
t/>(256(T + in the above inequality, we have 

suprj < 2V>(256(T + < 512^(512T) 

t<T 

for all sufficiently large T. This shows that 512-0(512£) is an upper rate 
function of Brownian motion on M under the probability P#i- The tech- 
nical point of passing from the average probability Pb 1 to the pointwise 
probability ¥ x is taken care of in the proof of the our main theorem below. 
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Theorem 4.1. Let M be a complete Riemannian manifold and let x 6 
M. Let B(R) be the geodesic ball on M of radius R and centered at z. Define 



and let tp be the inverse function of (p. There then exists a constant C such 
that Cip(Ct) is an upper rate function of Brownian motion X on M , that 
is, 



We have shown that Fb 1 (H) = 1. This shows that Cip(Ct) is an upper rate 
function for Brownian motion on M starting from the uniform distribution 
on the geodesic ball B\. Passing to a single starting point is easy. Let 



By the definition of the event If, it is clear that for any stopping time r, we 
have u) G H if and only if 9 t uj £ H; in other words, Lh ° T = Lh- It follows 
that h(z) = ¥ Z (H) = K z Lh is a harmonic function on M. On the other hand, 
we have < h < 1 and 



By the maximum principle for harmonic functions, we see that h must be 
identically equal to 1. □ 

The following special cases have all appeared in the literature (see the 
references cited in Section 1). They now follow from our main Theorem 4.1 
and all are now valid without any geometric restrictions. 

Corollary 4.2. Let M be a complete Riemannian manifold. Under the 
following volume growth conditions, ip is an upper rate function for Brown- 
ian motion on M : 




¥ x {d(X t ,x) < Ctp{Ct) for all sufficiently large t} = 1. 



Proof. Let 



H = {d(X t ,X ) < Cip(Ct) for all sufficiently large t}. 



h(z)=F z (H). 

Let 6 t : @>{M) ->■ @>{M) be the shift operator defined by 




(1) \B{r)\ < Cr D and ij){t) = CiVHiit; 

(2) \B(r)\ < e Cra (0 < a < 2) and *fj(t) = dt 1 ^ 2 -^ ; 

(3) \B(r)\ < e Cr2 and ip(t) = d exp(dt 2 lnt); 

(4) \B(r)\ < e Cr2lnr and ip(t) = exp(exp(GV)). 
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Proof. These upper rate functions follow directly from the main the- 
orem. Since the volume grows faster than the additional term lnlnr in the 
function <f>, these rate functions are the same as if the additional term were 
not there. □ 

Riemannian manifolds with slow volume growth are interesting test cases 
for our main result Theorem 4.1. Although, in general, slow volume growth 
corresponds to slow upper rate functions, our result will not yield upper 
rate functions better than y/i In In t once |-B(r)| < (lnr) 7 or some 7 > (see 
Remark 1.1). It also explains the condition v(r) > (lnr) 7 in the following 
result, due to Grigor'yan [2]. 

Corollary 4.3. Let M be a complete manifold such that \B(r)\ < v(r) 
for an increasing function v{r) > (lnr) 7 with some 7 > 0. Define R(t) by 

R(t) 2 

]nv(R(t)) ' 

Then CR{Ct) is an upper rate function for Brownian motion on M . In 
particular, if M has finite volume, then C y/tln In t is an upper rate function. 

Proof. With the lower bound for v(r), we have 

<f> f R )= f R rdr >c 1 f R rdr > ° 2R2 , 

Jq ]n\B(r)\ + lnlnr ~~ j 6 lnu(r) ~~ \nv(R) 

Therefore, the inverse function ip(t) < CsR{C$t). By Theorem 4.1, CR(Ct) 
is an upper rate function for some C. □ 

Remark 4.4. In all of the concrete cases we have mentioned thus far, 
upper rate functions are determined up to multiplicative constants. The 
question naturally arises as to whether we could have been more careful in 
our computations so as to recover the best constants in some cases, for in- 
stance, ip(t) = yj(2 + e)t mint for the standard one-dimensional Brownian 
motion. Such a precise upper rate function is impossible without further ge- 
ometric assumptions other than the volume growth. This can be explained 
by means of manifolds with power volume growth |-B(r)| < Cr D . According 
to Corollary 4.2(1), the corresponding rate function is ip(t) = C\\ft Int. By 
comparison with a Euclidean Brownian motion, we would expect a double 
logarithm instead of a single one. However, there are known examples show- 
ing that the above rate function with a single logarithm is indeed sharp up to 
a multiplicative constant (see Grigor'yan and Kelbert [5]). This is the reason 
why we have been somewhat cavalier about multiplicative constants in our 
proofs. It should be pointed out that these constants, denoted by C with or 
without subscripts, are universal; they do not depend on the manifold M 
(not even on its dimension). 
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